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ABSTRACT 

We present mid-infrared spectra of thirty two high redshift ultraluminous in- 
frared galaxies, selected via the stellar photospheric feature at rest-frame 1.6^m, 
and an observed-frame 24/im flux of >500/iJy. Nearly all the sample reside 
in a redshift range of {z) = 1.71 ± 0.15, and have rest-frame 1-1000/im lumi- 
nosities of 10^^'^ - lO^^'^L©. Most of the spectra exhibit prominent polycyclic 
aromatic hydrocarbon emission features, and weak silicate absorption, consis- 
tent with a starburst origin for the IR emission. Our selection method appears 
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to be a straightforward and efficient way of finding distant, IR-luminous, star- 
forming galaxies in narrow redshift ranges. There is however evidence that the 
mid-IR spectra of our sample differ systematically from those of local ULIRGs; 
our sample have comparable PAH equivalent widths but weaker apparent silicate 
absorption, and (possibly) enhanced PAH 6.2/im/7.7/im and 6.2/im/11.2/im flux 
ratios. Furthermore, the composite mid-IR spectrum of our sample is almost 
identical to that of local starbursts with IR luminosities of 10^'' — IO^^Lq rather 
than that of local ULIRGs. These differences are consistent with a reduced dust 
column, which can plausibly be obtained via some combination of (1) star for- 
mation that is extended over spatial scales of l-4Kpc, and (2) star formation in 
unusually gas-rich regions. 

Subject headings: infrared: galaxies — galaxies: active — galaxies: starburst — 
galaxies: evolution 



Introduction 



First discovered in the 1970s (IRieke fc Lowl Il972l ). Ultraluminous Infrared Galaxies 
(ULIRGs, those objects with l-lOOO/xm luminosities in excess of lO^^L©) are among the 
brightest objects in the low redshift Universe, with luminosities rivalling the bolometric 
output of quasars. Their immense luminosities, implying astrophysical processes oper- 
ating on extreme scales, have made ULIRGs the object of intense study over the last 
thirty years, and in the last decade or so a nascent consensus has developed on their 
nature. Local ULIRGs are probably 'composite' objects; most are powered mainly by a 
starburst, but a significant fraction (~ 45%) also contain an IR-luminous AGN. This is 



suggested from seyeral li nes of evidence, including opti cal/UV spectroscopv (IVeilleux et al. 



1999 



Farrah et a. 



Rigopoulou et al. 



199 



X-ray observations (IFranceschini et a] 



2005). mid- infrared sp e ctroscopy (|Genzel et a] 



Armus et al. 2007: Farrah et al. 



2003 



Ptak et a] 



sociated almost exc lusively with galaxy mergers (IFarrah et al.l l200ll : iBushouse et al 



1998 



Lutz et al 



2007bl : llmanishi et al.l 120071 ). and 



1998 



20031). Local ULIRGs are a 



1988 



Tacconi et al 



Veilleux et al.ll2002), and may be involved to some degree in triggeri ng QSOs (ISanders et al. 



2002; Kawakatu et al. 



so as- 



20021: 



2006 



ties of ULIRGs can be found in lSanders fc Mirabel! ( 
J2006h . 



Zauderer et al.ll2007f) . Reviews of the proper 



19961 ). and more recently in lLonsdale et al 



Initially, studies of ULIRGs were motivated mostly by their luminosities, as they are 
rare at low redshifts (with less than 50 examples known over the whole sky at z < 0.2), and so 
were not thought to play a fundamental role in the broader picture of galaxy evolution. This 
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perception changed abruptly however, when it was reahzed that ULIRGs are vastly more 
numerous at high redshift. First hinted at by spectroscopic foUowup o f Infrared Astronomica l 



Satellite (IRAS) surveys (iHacking et al.l 119871 : iLonsdale et al.l 119901 : ISaunders et al.l Il990l ) 



which showed strong evolution in the ULIRG luminosity fun ction with redshift, this was 



confirmed by surveys with the Infrar ed Space Observatory fI SQ iRowan- Robinson et al 



Dole et al 



2000 



2001 



Scott et alll2002l : iBorvs et al. 



Verma et al.l l2005h and in the sub-mm (Hughes et al 



1998 



1997 



Bales et al. 



20031 : iMortier et al.ll2005l ). which showed that there were 



several hu ndred ULIRGs per square degree at z > 1. Later surveys with the Spitzer space 
telescope (jWerner et al.ll2004j ) showed that ULIRGs contributed an increasing fraction of the 
comoving bolometric energy density with increasing redsh ift, reaching a level co mparable to 
the contribution from lower luminosity systems at 2; > 1 (ILe Floc'h et al.ll2005l . their figure 
14). Distant ULIRGs seem to be similar to their low red shift counterparts, i n that they 



appear to be powered by both starburst and AGN activitv (IFarrah et al 



2OOII2OO4J : [Alexander et al. 



mostly mergers (IFarrah et al. 



2005 



2002a 



Takata et al. 2006: Valiante et al. 



Chapman et al. 



2003 



Bridge et al. 



2002bl : ISmail et al. 



2007), and are probably 



20071 ). 



The much higher space density of ULIRGs in the high redshift Universe compared 
to locally implies that they play an increasingly important role in galaxy evolution with 
increasing redshift; for example, their high star formation rates make distant ULIRGs ex- 



cellent candidate s for being the rapid growth phases of massive 



( IScott et all I2OO2I : iRocca-Volmerange et al.l l2004l : ISwinbank et al. 



> 2L *) elliptical galaxies 



20061 ). The sheer number 



of ULIRGs at high redshift however also raises some serious problems, notably for theories 
for the formation of large-scale structures. These theories usually invoke 'biased' hierarchi- 
cal buildup, in which overdense halos in the dark matter distribution undergo successive 
mergers to build haloes of increasing mass, with galaxies forming from t he baryonic matter 
in these haloes (ICole et alll2000l : iGranato et all l2000l : iHatton et alll2003l ). In this scenario, 
ULIRGs, being the likely antecedents of massive systems, should be found i n overdense peaks 



in the underlying D M distribution, and observations seem to bear this out (IBlain et al. 



2004 



Farrah et al.H2006al jbll Despite this, the prevalence of ULIRGs at high redshift, and their 
very high inferred star formation rates, cause major problems for models as they find it 
difficult to concentrate enough baryons in suffic iently small volumes, an d a diverse variety of 



Bauffh et al 


2005: 


Bower et al. 


2006) 



There is therefore a pressing need for both an efficient census of ULIRGs at 2 > 1, 
and a thorough understanding of how these distant ULIRGs form stars. Spitzer is ideal 
for undertaking such studies; the Infrared Array Camera (IRAC, Fazio et al.l l2004j ) and 



Multiband Imaging Photometer for S pitzer (MIPS. iRieke et al.ll2004j ) imaging instruments 



and the Infrared Spectrograph (IRS. iHouck et alll2004j ) all offer dramatic improvements in 



-4- 



sensitivity and resolution compared to previous generation facilities. In this paper, we use the 
IRS to study 32 ULIRGs at 2 > 1, selected using IRAC and MIPS photometry. We quantify 
the efficiency of the selection method and use mid-IR spectral features to set constraints on 
the nature of their star formation. An overview of the sample selection and a description of 
the observations is given in ^ Results are summarized in ^with discussion and conclusions 
given in §l]and ^respectively. We assume a spatially flat cosmology with Hq = 70 km s~^ 
Mpc~^, n = 1, and flm = 0.3. Unless otherwise stated, the term 'Infrared (IR) luminosity' 
refers to the total luminosity integrated over 1-1000/im in the rest-frame. 



2. Analysis 



2.1. Sample Selection 



The basis of our selection is an apparent emission feature at rest-frame 1.6/im, produced 
by thermal emission from late-type stars and enhanced by an apparent emission feature 
caused by H" ions in the atmospheres of giant star^. T his 1.6/im 'bump' is almost ubiquitous 
in th e integrated spectra of galaxies. Previous authors (ISimpson fc Eisenhardtlll999l : ISawicki 
2OO2I ) have investigated the use of this feature in deriving photometric redshifts for galaxies 
using the Spitzer IRAC channels, and conclude that it is likely to be effective. A 1.6/im 
bump also implies that any AGN continuum is not signiflcant in the rest-frame near-IR, 
although it does not rule out the possibility of a low luminosity and/or very obscured AGN. 
Therefore, combining the presence of the 1.6/im feature with a prior on IR luminosity should 
select high redshift ULIRGs that harbor active star formation. A complete description of the 



selection method can be found in Lonsdale et al (2008, in preparation). Angn 
measur ements for large samples of ' bump' selected ULIRGs can be found in 



fl2006ar ) (see also lOliver et al.ll2004l : IWaddington et aL 



2007 



ar clustering 



Farrah et al. 



de La Torre et al.l 120071 ) . and 



millimeter- wave fluxes and analysis can be found in Lonsdale et al. ~ (|2007b). Other studies 



that include small numbers of 'b ump' sources inc lude IRS spectroscopy (IWeedman et al. 



20061 ) ■ and optical spectroscopy (IBerta et al.l 120071). O ther studies of the IRAC colors of 
SWIRE galaxies can be found in bavoodi et all J2006al jbll 



We select our sources from the SWIRE (ILonsdale et al.l |2003| ) Lockman fleld, which 
spans 10.2 square degrees, and has been imaged in all seven IRAC and MIPS bands to 5a 
depths of 4.2/iJy, 7.5/iJy, 46.1/tJy, 47.3/iJy, 210/iJy, 18mJy and lOSmJy respectively. We 



is a m ajor source of opacity in the atmospheres of such stars, but shows a minimum in opacity at 
around 1.6/im (John IQSSh . leading to a broad apparent emission feature at this wavelength 
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first excluded all sources brighter than nir = 23 to remove local and low redshift galaxies, 
and all sources fainter than 500/iJy at 24/im to ensure that spectra could be obtained with 
the IRS in reasonable integration time^. We then selected all sources that exhibited excess 
emission in the 4.5/xm channel, i.e. where /s.g < /4.5 > /s.s and f^ ^ > fs.o (or 'bump 2' 
sources). This resulted in a total of 350 sources. We then randomly selected 31 objects 
for observation with IRS. We also included one further object in the SWIRE ELAIS Nl 
field that satisfies all of our selection criteria, and was observed as part of a complementary 
program. 

The complete sample of 32 objects is listed in Table [H The IRAC and MIPS fluxes in 
this table are from the latest SWIRE data release (the 'SWIRE3' catalogs). The sources were 
however originally selected using the previous generation, 'SWIRE2' catalog. The changes 
to the IRAC pipeline between these two catalog versions resulted in flux differences between 
the catalogs of only a few percent, but our original selection did not demand that the bump 
was detected with any statistical robustness, and in some cases the band 2 excess was not 
significant within the errors. As a result, some of our sample no longer display an excess in 
the 4.5/xm channel using the SWIRE3 catalogs. The IRAC fluxes of all of the sample are 
however consistent within the errors with an excess in the 4.5//m channel, and for simplicity 
we treat all of the sources in Tabled] as 'bump 2' sources for the remainder of this paper. 



2.2. Observations 

Each source was observed with the Long-Low module onboard the IRS (10.5" x 168", 
5.1" pix-\ R~ 57 - 126) in both first (LLl, 19.5//m - 38.0//m) and second (LL2, 14.0^m 
- 21.3/im) order, as part of Spitzer program 30364 and (for source 32) program 15. The 
targets were placed in the center of each slit using the blue peak-up array. Each target was 
observed for eight ramp cycles, with an individual ramp time of 120 second^. 

The data were processed through the Spitzer Science Center's pipeline reduction software 
(version 15.3), which performs standard reduction tasks such as ramp fitting and dark current 
subtraction, and produces Basic Calibrated Data (BCD) frames. Starting with these frames. 



^The implie d 24/optical color-cut is substantially more lenient than that of Houck et al. 20051 and 



Yan et al.ll2007l which would have excluded an Arp 220 type galaxy at z=2. 

^Details for each observation can be found by referencing the AOR keys given in Table [1] within the 
Leopard software, available from the Spitzer Science Center 
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we flagged rogue and otherwise 'bad' pixels using the zrsc/earfl tool, which uses a pixel mask 
for each campaign to first flag and then replace rogue pixels. The individual frames at each 
nod position were t hen combined into a single image using the SMART software package 
(iHigdon et al.l 120041 ) . Sky background was removed from each image by subtracting the 
image for the same object taken with the adjacent order (i.e. 'order-order' sky subtraction). 
One-dimensional spectra were then extracted from the images using the SPICE software 
package using 'optimal' extraction and default parameters for the Long-Low modules. This 
results in separate spectra for each nod and for each order, all of which must be combined 
to give the final spectrum. The spectra for each nod were inspected; features present in only 
one nod were treated as artifacts and removed manually. The two nod positions were then 
combined. In most cases the two orders differed in flux by ~ 20%; we first scaled LL2 to 
LLl, and then scaled the resulting spectrum to the 24//m flux. Finally, the first and last 
4 pixels on the edge of each order, corresponding to regions of decreased sensitivity on the 
array, were removed to give the final spectrum for each object. 



3. Results 

The spectra are presented in Figures [1] and [2l Redshift data are presented in Table |2l 
and measurements of spectral features are given in Table [31 

We derive spectroscopic redshifts from broad emission features at 6.2/im, 7.7/im, 8.6//m, 
11.2/im and 12.7//m, attributed to bending and stretching modes in neutral and ionized PAH 
molecules (the 12.7/im feature likely also contains a contribution from the [NeII]A12.814 fine- 
structure line), and/or a broad absorption feature at 9.7/xm arising from large silicate dust 
grains. In most cases (29/32) an unambiguous redshift can be determined from the PAH 
features; the associated uncertainty on these redshifts is governed by the variations in PAH 
peak wavelengths seen in local galaxies, and is of order Az = 0.02. In three cases however 
(sources 11, 15 and 27) the spectral features cannot be reliably identified. For these three 
sources we derive tentative redshifts based on what are plausibly PAH or silicate features. 
In these cases the error on the redshift is larger, Az ~ 0.3, and the redshift itself should 
be treated with caution. A histogram of the redshifts is plotted in the left panel of Figure 
[3l All of the sample lie in the redshift range 0.8 < z < 2.3, with all but three lying within 
1.4 < z < 1.9. The redshift distribution is well fitted by a Gaussian centered at z=l.ll and 
with a Full Width at Half Maximum of Az~0.35. 

An accurate measure of the IR luminosities is difficult due to the paucity of data at far- 



''This tool is available from the SSC website: 



http : / / SSC .spitzer . caltech.edu 
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IR wavelengths. Therefore, we initially simply estimate the range in IR luminosity within 
which our sample likely resides, and defer estimates of IR luminosities based on the IRS 
spectra to ^4.3[ To do this, we no rmalize a suite of radiative traiisfer m odels for starbursts 
fafstathiou et aliboooh and AGN jEfstathiou fc Rowan- RobinsonlEoQsl ) to the 2Aijm fluxes, 
and reject those models that exceed the IRAC fluxes, or 70;um or 160/im upper hmits. The 
IR luminosities derived from the remaining models then constitute the probable range in IR 
luminosity of our sample. The resulting IR luminosity range is 10^^'^ — placing our 

sample at the upper end of the IR luminosity range seen for low-redshift ULIRGs 

The sample is faint at 24/im, hence our spectra are of low signal-to-noise, and detailed 
spectral diagnostics are not possible. Most of the spectra however show one or more PAH 
features at 6.2/im 7.7/im 8.6/im, 11.2/im and (depending on redshift) 12.7/im. The PAH 
fluxes and equivalent widths (EWs) were computed by integrating the flux above a sp line 
interpolated local continuum fit (for a description of the method see lBrandl et al.ll2006l and 
Spoon et al.l 120071 ). Due to the low S/N and restricted wavelength range of the spectra, we 
cannot correct for water ice and/or aliphatic hydrocarbon absorption, however the effect of 
this lack of correction is likely to be insignificant. The errors on the EWs are large because 
the method used to measure them requires a spline fit to the underlying continuum, and the 
continua of our sample are in most cases only weakly detected. This method does however 
offer the advantage that our E Ws can be compared directly to those o f local IR-lurninous 



Armus et al. 


2007; 


Spoon et al. 


2007) 



We derive star formation rates from the cornbined lu minosity of the 6.2/im and 11.2/im 
PAH features, via the formula from iFarrah et al.l (j2007bl ) : 



M^= 1.18 X 10-^^Lf 



(1) 



where Mq is the star formation rate in solar masses per year, and Lp is the PAH 6.2/im + 
11.2/im lum inosity in Watts. Th e resulting star formation rates are listed in Table |3l As 
described in lFarrah et al.l (j2007bl ). the errors on these star formation rates are of order 50%. 



We measured the strengths of the silicate features via: 



S. 



sil 



Fcont{Q-7fim) J 



(2) 



where Fobs is the observed flux density at rest-frame 9.7/tm, and Fcont is the underlying 
continuum flux at rest-frame 9.7/im deduced from a spline fit to the continuum on either 
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side. A complete description o f the method used t o mea sure Ssu can be found in lSpoon et al. 
(120071 ) ; iLevenson et al.l (120071 ) and ISirocky et al.l (120071 ) . The sihcate feature is not detected 
(at > 3a) in any object individually. We also see no evidence for molecular Hydrogen lines 
or any fine-structure emission lines in the individual spectra. 

The average spectrum for all 32 objects is shown in Figure HI All of the PAH features 
between 6.2/xm and 12.7//m are prominent. There is also evidence for two fine-structure lines; 
the 12.7^m PAH feature appears asymmetric, likely due to the presence of [NeII]A12.814, and 
there is a weak d etection of what may be [ArII]A6.99. B oth these lines are commonly seen 
in local ULIRGs (iGenzel et allll998l : iFarrah et al.ll2007bl ). Conversely, the average spectrum 
shows no evidence for significant apparent silicate absorption. We also see no signs of the 
H2,S(3)A9.662 and H2,S(2)A12.275 lines, which are both common, though weak, in local 
ULIRGs, or the [SIV]A10.511 line. These results do not change appreciably if we selectively 
coadd the spectra, examples of which are also shown in Figure HI The average of the ten 
spectra with the largest PAH 6.2yum EWs shows no apparent silicate absorption, no molecular 
Hydrogen features, but strong PAH features with an asymmetric 12.7/im PAH feature. The 
average of the ten spectra with the 'strongest' silicate absorption features is similar; there is 
a hint of absorption, but it is not detected with any statistical significance. 



4. Discussion 



4.1. Redshifts and luminosities 



The redshifts of our sample have a remarkably tight distribution. All of the sample lie 
in the range 0.8 < z < 2.3, with all but three lying within 1.4 < 2; < 1.9. If we exclude 
these three outliers, then the mean redshift for our sample is (z) = 1.71 ± 0.15 (assuming 
FWHM= {2y/2h}2)a). This is significantly tighte r than the redshift d istribution of sub- 
mm selected galaxies, which have (z) = 2.2 ± 0.5 ( Chapman et al. 20051), and the redshif t 
distribution for Lyman break galaxies, which have {z) = 2.96 ± 0.29 (ISteidel et al.ll2003l ). 
though a larger sample would be needed to confrm this. Our redshift dis tribution is muc h 
tighter than that of sources found by selecting on 24//m/0.7^m flux ratio (lYan et al.l 120071 ). 
which span 0.5 < z < 3.0. Combined with the constraint that the IR luminosities of our 
sample must (for any reasonable SED shape) be in excess of 10^^-^Lq, we conclude that our 
selection method is an excellent way of finding distant ULIRGs in specific, narrow redshift 
intervals. This narrow redshift range likely arises due to a combination of two factors; the 
demand of an excess flux at observed frame 4.5/Ltm favors sources with a rest-frame 1.6//m 
opacity feature at 1.4 < z < 2.0, and the demand that sources be bright at observed-frame 
24//m biases towards sources whose redshifts place the 7.7/im and 8.6//m PAH features 
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within the MIPS 24/im band. This method offers many advantages over other methods for 
finding high redshift ULIRGs. For example, it avoids all of the observational difficulties in 



assigning optical counterparts to sub-mm selected ULIRGs (see for example iDunlop et al. 



(120041 ) )■ Remarkably, we have arrived at this redshift distribution without demanding that 
the 'bump' be detected with any statistical robustness, and it is notable that none of the 
seven sources that 'lost' their bump in the 4.5/im channel due to reprocessing (Table [2]) are 
unusual in redshift or spectral shape. 



Given the difficulties in obtaining accurat e redshifts for distant ULIRGs (ISimpson et al. 



2004J : IChapman et al.l l2005l : iBerta et al.l 120071 : lAretxaga et al.l 120071). we examin e the accu- 
racy of photometric redshifts for our sample obtained by iRowan-Robinson et al.l (120071 ) and 
listed in Table [21 Given the high redshifts and limited optical photometry, the photomet- 
ric redshifts are surprisingly good. 'Rehable' solutions (those with Xred < ^0 detected 
in at least three of the Ugriz bands) are obtained for 14/32 objects; these are plotted 
in the right panel of Figure [31 Of those 14, four are classified as 'catastrophic' failures, 
i.e. those where the photometric redshift deviates from the spectroscopic redshift by more 
than 6% in log(H-z) (although two of these four are sources 11 and 15, which have un- 
certain spectroscopic redshift s), while ten are classified as successes. We conclude that the 
Rowan-Robinson et al.l (120071 ) photometric redshift method remains effective even for faint, 
high redshift sources. Finally, we note that, of the seven objects in Table [2] with ambiguous 
bump classifications, six have photometric and spectroscopic redshifts that agree very well. 
It seems therefore that, for the sources with an ambiguous bump classification, most are 
indeed 'bump 2' sources. 



4.2. What powers the infrared emission? 



The IRS spectra provide the best diagnostics of the power source behind the IR emission. 
The presence of prominent PAH features in 29/32 objects is strong evidence that > 90% of 
our sample harbor powerful starbursts. The star formation rates computed using Equation 
[H and listed in Table [3] are all extremely hi gh. If we extrapo late these star formation rates 
to IR luminosities via the relation given by Kennicutt ( 1998 1^. then these luminosities lie in 
the same range as those determined in §31 We infer that dust heated by star formation can 
account for all of the IR emission from the sources that show strong PAHs, and therefore that 
most of our sample are consistent with being pure starbursts. Interestingly, the strong PAHs 



KennicuttI (jl998[ ) derive a relation between star formation rate and rest-frame 8-1000/im luminosity 



rather than 1-1000/xm luminosity, but the difference in wavelength range is unlikely to be significant 
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and high degree of homogeneity of our sample make them similar to sub-mm selected ULIRGs 
( iMenendez-Delmestre et al.ll2007l : lPope et alll2007l ). suggesting that we may be drawing from 
the same parent population. 

We cannot, however, rule out the presence of an AGN in any of our sample. It is 
unlikely that an AGN contributes significa ntly in the rest-frame near-IR, as we see no powe r 



law continua in the IRAC bands (see e.g. lAlonso-Herrero et al.ll2006l : iPoUetta et al.ll2006l ) 



but it is possible that some of our sample contain mid-IR luminous AGN, even in those 
objects with prominent PAHs. An example of such an object at low redshift is NGC 6240, 
which shows st rong PAHs accomp anied by an AGN continuum and the [NeV]A14.322 fine- 
structure lin^ ( Armus et al. 2006 ). Furthermore, previous authors have noted correlations 
between shal l ow (S.^u < 0.8) or deep {Sgu > 2.4) silicate features, and the presence of AGN; 
Farrah et al.l ( 2007bl) show th a t ULI RGs with shallow silicate features usually also show the 



[NeV]A14.322 line, iHao et al.l (120071) find tha t QSOs and Seyfert 1 galaxies often show the 
silicate feature in emission, and ISajina et al.l (j2007bl ) find that many high redshift ULIRGs 
with deep silicate absorption features also show powerful radio jets. The shallow silicate 
features of our sample are therefore consistent with (though not supportive of) the presence 
of an AGN. 

Finally, we examine the IRAC colors as a diagnostic of the IR emission. In Figure [5] we 
plot 5.8/im /3.6//m color again st 8.0/um/4.5/im color, together with the selection 'wedge' for 
AGN from iLacy et al.l (120041 ). defined empirically using the Spitzer colors of SDSS QSOs. 
This diagnostic is not, on its own, particularly useful, as our 4.5/im 'bump' selection already 
excludes much of the AGN selection region, but it serves as a useful check. Nearly all of 
our sample lie outside the AGN selection region, and are consistent with the positions of the 
bulk of the galaxies in the main Spitzer First Look Survey (FLS) dataset. Our sam ple lies 
some distance from the SDSS QSOs and Syl galaxies in figure 1 of iLacy et al.l (120041 ). Based 
solely on the IRAC colors, the rest-frame near-IR emission from our sample is consistent 
with arising from starlight. 



^[NeV] A14.32 has Eion — 97.1 eV, and models indicate that it is unHkely to be strong in galaxies without 
an AGN (jVoitl 119921 ) . Qbservationally , [NeV]A14.32 is occasionally seen in planetary nebulae and young 
supernova remnants (jOliva et alJll999l ) , but is not seen in the integrated spectra of starburst galaxies. It is 
however strong in the spectra of AGN ( Lutz et aJ 1998h 



4.3. The nature of star formation in distant ULIRGs 



4.3.1. PAH flux ratios 

The prevalence of starbursts in our sample means we can investigate the nature of star 
formation in high redshift ULIRGs. We start by considering diagnostics based solely on 
the PAH features. In Figure [6] we plot PAH luminosities against PAH flux ratios. The 
6.2/im and 7.7/im PAH features are extraordinarily luminous, lending further weight to the 
conclusion that the starbursts in our sample are more intense than seen in any local [z < 0.2) 
ULIRG. The dispersion in PAH luminosity is perhaps slightly greater for our sample, but 
not signiflcantly so given the sizes of the error bars. 

We see a more interesting picture if we consider PAH color-color plots (Figure [7j). Con- 
sidering flrst only the low redshift sources; local IR-luminous galaxies show a distinct 'en- 
vironmental segregation' with IR luminosity, corresponding to elevated 6.2/im/11.2/im and 
7.7^m/11.2/im PAH ratios in ULIRGs compared to lower luminosity starbursts and AGN. 
Such a segregation is als o seen in HII regions and Photo-Dissociation regions, and the rea- 



sons for it are not clear (jVermeij et al.l 120021 : iLebouteiller et al.l 120071 ) . The two most likely 



reasons are PAH ionization and PAH molecular structure. A higher fraction of ionized 
PAHs gives rise to stronger PAH features at shorter wavelengths, as does a higher fraction of 
open, uneven PAH molecules compared to smooth, compact PAH molecules (see figure 9 of 



Hony et al.l (120011 ) for diagrams of these PAH features). The drivers behind PAH molecular 
structure are unknown, but PAH ionization fractions can be affected in two ways; (1) a harder 
interstellar radiation field (ISRF) would ionize more pahA, and (2) decreasing the electron 
density would decrease the recombination rate, leading to a higher fraction of ionized PAHs. 
It is therefore likely that local ULIRGs have either a harder ISRF, a lower electron density, 
or a higher fraction of open, uneven PAH molecules than do lower luminosity systems. 

Turning now to the bump sources; our sample appears to lie preferentially toward one 
'edge' of the local ULIRG distribution, corresponding to even further elevated 6.2/im PAH 
emission compared to the 7.7/im and 11.2/im features, though the error bars are large even 
for the average bump point. These two diagrams alone suggest that the star formation in our 
sample may differ in some way from that seen in local ULIRGs, although it is difficult to say 
what the physical origin of those differences could be. One possibility is a lower extinction 
level in our sample compared to the mean for low redshift ULIRGs. 



^though ISmith et al.l (|2007r ) note that a hard ISRF from a low-luminosity AGN may actually suppress 
PAHs in the 5-8/im range 
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4-3.2. PAH fluxes & Silicate strengths 



We move on to consider diagnostics based on both the PAH emission features and the 
9. 7jj,m sihcate a bsorp tion feature. In Figure [8] we plot our sample on the 'fork' diagram 
of ISpoon et al.l (120071 ). Our sources occupy a small region on the lower right of this plot, 
corresponding almost entirely to class IC, with a small number in classes 2C and 2B, though 
the large error bars on the PAH EWs mean our sources are in many cases also consistent 
with class lA/lB. These classes label our sources as PAH dominated, with relatively weak 
apparent silicate absorption. To check whether this is just a chance similarity between these 
specific spectral features, we plot in Figure M the averaged spectrum from Figure H] against 
the averaged spectra for the classes in the fork diagram. Clearly, our averaged spectrum 
closely resembles the average spectrum for the class IC sources. 

Figure [8] raises interesting questions on the nature of starburst activity in our sample. 
As described in §2.H our selection method biases towards ULIRGs harboring ongoing star 
formation, and selects against those ULIRGs containing a rest-frame near-IR luminous AGN. 
Our sample should therefore resemble the majority of local ULIRGs, a s most low redshift 
ULIRG are powered mainly by star formation (e.g. iFarrah et al.l 120031 ). Furthermore, our 
demand for an excess flux in one of the IRAC channels should not remove many 'local type' 
ULIRGs from the sample, as many local ULIRGs have a rest-frame SED that either has a 
near-IR excess, or is not obviously a power law. Moreover, previous studies of high redshift, 
L/ij ~ IO^^Lq ULIRGs selec ted using Sp i tzer da ta have shown that they seem to resemble 
local ULIRGs; for example, ISajina et al.l (l2007al ) find that z > 1 ULIRGs selected on the 
basis of red 24/im/8/im and 24/im/0.7/im color often show PAHs accompanied by an AGN 
continuum and in some (~ 25%) cases strong silicate absorption, and their figure 13 appears 
to place the bulk of their sample in the upper branch of Figure Ei We might therefore expect 
that our sources should be distributed mainly along the upper branch in Figure [H along 
with the majority of local ULIRGs, and previous small samples of Spitzer selected high 
redshift ULIRGs. This however is not the case; our sources are located in a specific region, 
classifying them as starbursts with (relatively) weak apparent silicate absorption, and almost 
completely avoid the upper branch in Figure [8], where most local ULIRGs reside. 

There is no a priori reason why our relatively crude selection method should select 
against ULIRGs with strong silicate absorption, and the almost complete absence of such 
absorption in our sample is peculiar, for two reasons. First, many local ULIRGs have a 



near- IR excess coupled with strong silicate absorption (ISpinoglio et al.l Il995l : lArmus et al. 



20071 ). Second, the 24/im fiux cut (in combination with the demand for a rest-frame l.G^um 
'bump') selects in favor of sources with strong PAHs, but among loca l ULIR Gs, systems with 
strong PAHs but weak silicate absorption are rare - e.g. iDesai et al.l (120071 ) examine ULIRG 
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local spectra medianed by optical spectra type, and find strong PAHs and prominent silicate 
absorption in the HII and LINER composites, and weak PAHs and weak silicate absorption 
in the Seyfert composite. 

If our sample do not resemble local ULIRGs, then what do they look like? We com- 
pare in Figure [TOl o ur average spectrum ^ yith three different starburst templates; Arp 220 
(lArmus et al.l 120071 ). the nucleus of M82 (ISturm et al.l 120001 ). and the average spectrum of 
the local starburst galaxies from iBrandl et al.l |2006| (which have IR luminosities between 



IO^^Lq and IO^^'^Lq). Clearly, the Brandl et al spectrum is the closest match to our 
spectrum, even though the parent samples differ in IR luminosity by over two orders of 
magnitude. Both show strong PAHs with an asymmetric 12.7/im PAH profile, and negli- 
gible silicate absorption. Interestingly, another class of high redshift gala xy that tends to 
show strong PAHs but weak sili cate absorption are the sub-mm galaxies (jPope et all 120071 : 
Menendez-Delmestre et al.l 120071 ). The M82 nuclear spectrum is also a good match, except 
for slightly stronger silicate absorption. The Arp 220 spectrum on the other hand is a poor 
match, with much stronger silicate absorption. 

The similarity of our composite spectrum to that of local starburst galaxies allows us 
to estimate the IR lumino sities for our sampl e using the relation between IR luminosity and 
6.2/im PAH luminosity in lBrandl et al.l 12006 : 



log{Lin) = 1.13 x log {F^.2DI) (3) 

where Fq,2 is the flux of the PAH 6.2^m feature in units of 10~^^ W cm~^, Dl is the luminosity 
distance in Kpc, and Lm is in units of solar luminosities. The IR luminosities for our sample 
computed using Equation [3] are presented in Table |31 and are in good agreement with the 
IR luminosities estimated in ^ 

We explore the robustness of the conclusion that our sample have systematically weaker 
apparent silicate absorption than is seen in local ULIRGs using a different approach; fitting 
their br oad band optical through IR S EDs (i.e. the U griz, IRAC and MIPS data) with the 



code of iRowan- Robinson et al.l (l2007l ). with the redshifts fixed at the spectroscopic values, 
and comparing the results from these fits with the IRS spectra. We consider only those 
objects (4, 12, 13, 15, 17, 21, 29, 31) with detections in at least six optical through MIPS 
bands, and specifically do not fit to the IRS spectra, as we wish to test the conclusions 
independently. The limited available photometry (less than nine points in all cases) renders 
this approach crude compared to the analysis of the IRS spectra, but it serves as a useful 
check. The results from this fitting are shown in Figure [TTl In four cases (4, 13, 21 29), 
an M82 template gives an excellent fit to the photometry, and is well matched to the IRS 
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spectrum. In two cases (12 and 17) an M82 template gives a good fit to the photometry, 
but predicts a higher fiux at rest-frame 6/xm than is seen in the IRS spectrum, which may 
indicate a steeper mid-IR continuum. In the final two cases (15 and 31) there is marginal 
evidence for strong silicate absorption, consistent with an Arp 220 template, though the low 
S/N of the IRS spectra means we cannot rule out an M82 template either. Overall, the 
results from the template fitting agree with the results from the individual IRS spectra, and 



from the comparison between our average spectrum and Arp 220, M82, and the lBrandl et al. 
20061 composite in Figure dUl 



All of the evidence points towards an unexpected result; despite the fact that our sample 
are ultraluminous in the IR, their mid-IR spectra have systematically lower silicate absorp- 
tion than low-redshift ULIRGs, and instead resemble the mid-IR spectra of local starbursts 
that are over two orders of magnitude less luminous in the IR. What makes this even more 
puzzling is that high redshift ULIRGs appear to be physically similar to local ULIRGs in 



many ways 



: t 



Bridge et al. 



ley sh ow evidence for b eing mergers ( Farrah et al.ll2002al : IChapman et al.ll2003 



flTacconi et al.ll2006l : iGenzel et al.ll2006l ). 



20071. though see also iGenzel et al.l l2006l ) between massive, gas rich systems 



We can think of six possible reasons to explain this difference: 

1 - IR Luminosity dependent silicate strength: Our sample is significantly more IR-luminous 
than most local ULIRGs, so we would see weak apparent silicate absorption if compact, dusty 
starbursts and AGN exhibit weaker silicate absorption as they increase in IR luminosity, but 
remain otherwise identical. While we cannot rule this possibility out, we can think of no 
plausible physical reason why this should be the case. 

2 - A specific evolutionary stage: We may be selecting high redshift ULIRGs at a specific 
evolutionary stage, when a compact but relatively old starburst (> 20Myr if we assume an 
exponentially decaying burst) becomes less obscured as supernova winds clear dust from the 
star-forming regions. While this is in principle plausible, our selection method is relatively 
crude compared to others cited as selecting distant ULIRGs at specifi c evolutionary stages , 
which rely on the presence of specific fine-structure absorption lines (IFarrah et al.ll2007al ). 
It seems unlikely that our method should find all 32 of our sample at this specific stage. 

3 - Silicate poor dust: Weak apparent silicate absorption could arise from dust that is unusu- 
ally poor in silicates. If the dust is made mainly in Asymptotic Giant Branch (AGB) stars 
( lGehrzlll989l ). then silicate poor dust could be made in two ways; a low metallicity (Z< 0.04) 
would suppress the formation of silicate dust, and a steep or 'top-heavy' stellar IMF may 
shift more AGB stars in to the carbon dust producing mass range , thereby reducing th e pro- 
duction of sihcate dust (jBlanco et al. 1978 : Renzini fc Voli 1981; Lagadec et al. 2007). The 
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meta llicities of distant ULIRGs are not well constrained (IShapley et al.l |2004J : iLiang et al. 
20041 ). but it seems unlikely that the metallicities of our sample could be lowered to the 



extent required, as low metallicit ies also suppress PAH features (lEngelbracht et al.l 12005 



Smith et al.l 120071 : IWu et al.l 120071 ) . A top heavy IMF is more plausible, an d is consistent 



with the demands from structure formation models (e.g. iBaugh et alll2005l ). Interestingly, 
if the IMF were top heavy, it would result in a harsher radiation field, which (as described 
in §4.3. ip may result in enhanced 6.2/im PAH emission compared to the longer wavelength 
PAH features, a result hinted at in Figures 171 and [TUl There is however no direct evidence for 
a top-heavy IMF in distant ULIRGs. Moreover, some authors have suggeste d there is no evi- 
dence for evolution in dust properties with redshift, at least up to z ~ 1 
2OO5I : ISavouzet et all[2007l ). 



[e.g. 



Clements et al. 



4 - An intrinsically weak underlying continuum: If the underlying near/mid-IR continua in 
our sample are weak, then we would see three effects; the rest-frame 1.6/zm feature would 
be prominent, the PAH EWs would be higher, and we would see weak silicate absorption as 
there would be no background continuum upon which to make such absorption apparent. A 
weak continuum would arise if our sample harbored lower luminosity AGN than those seen in 
local ULIRGs. In this sc enario the ab s ence of silicate absorption, and the differences between 
our sample and that of ISajina et al.l (l2007al ). arises due to selection bias, as our selection 
method favors systems harborin g active star formatio n. It is also possible that the differences 
between our sample and that of ISajina et al.l (l2007al ) arise because our sample is fainter at 
observed-franie 24)um (most of our sample fall below the 24/im flux cut of 0.9mJy used by 
Sajina et al.l (l2007al )) which, if true, would imply a 'transition' in the mid-IR properties of 
high redshift ULIRGs at /24 < ImJy. This possibility is simple, but does suffer from three 
drawbacks. First, as described in §2.H we do not demand that the 1.6/im bump is detected 
with any statistical significance, and many of our sources have a weak 1.6^m bump. Our 
selection method therefore does not demand that the underlying continuum in the rest-frame 
near-IR is negligible, and therefore by extension does not demand that the rest-frame mid-IR 
continuum is absent either. Second, as can be seen in Figures [1] and [2], several of our sample 
have significant mid-IR continua. Third, several local ULIRGs with strong mid-IR continua, 
in some ca ses accompanied by si licate absorption, also show a near-IR 'bump' (e.g. Arp 220, 
NGC6240, ISpinoglio et al.lll995l ). hence the presence of a bump does not necessarily select 
against objects in the upper branch of Figure [81 

5 - Gas-rich mergers: A simple way to reduce the apparent silicate strength is to lower the 
total extinction, which can be achieved by reducing the amount of dust in the starburst 
regions, or equivalently by increasing the amount of gas so that the same level of star 
formation can be sustained for a smaller dust mass. This scenario is compatible with the 
presence of strong PAH emission, as changing the dust-to-gas ratio mainly affects the total 
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dust column to the central ionizing source, whereas PAH emission is thought to originate 
from the surface layers of star-forming regions. There is some support for this possibility 
from millimeter interferometry of sub-mm selected ULIRGs at h igh redshift, which shows 
that they may be unusually gas-rich compared to local examples ( iTacconi et al.ll2006l ). and 
also from models for the star formation history of the Universe, which suggest a lower mean 
Ay at z ~ 1.5 compared to locally (e.g. iRowan-RobinsonI |2003| . their figure 20). 



6 - Starburst geometry: The final possibility is a different star-forming region geometry, such 
that the total dust column is reduced. This can be achieved in two ways; either make the 
star formation diffuse and extended over several Kpc instead of the sub-Kpc scale starbursts 
seen in local ULIRG^, or by distributing the star formation in multiple compact dusty 
star forming regions spread over a few Kpc, instead of a single sub-Kpc nuclear starburst. 
There is observational evidence for both these 'modes' of star formation in distant starburst 
galaxies; there is evidence for extended star formation from hi gh resolution radio observations 
( IChapman et al.ll2004j : iMuxlow et al.ll2005l ) and Lya imaging (IGeach et al.ll2005l ) , while other 
authors have found evidence for star formation in multi ple small regions di stributed across 
a few Kpc in forra ing massive galaxies at high redshift (IGenzel et al.ll2006l ). Furthermore, 
Daddi et al.l (120071 ) have recently suggested that disklike, 'inefficient' star formation may be 
prevalent among distant massive galaxies. There is however an observational limit on the 
spatial extents of these starbursts; high resolution millimeter imaging of small samples of 
sub-mm selected ULIRGs at high redshift hav e shown that the st arbursts in these systems 
cannot be much more than 4Kpc in diameter (ITacconi et al.ll2006l ). 



It is difficult to discriminate reliably between these six possibilities, as the IRS spectra 
cannot formally rule any of them out. The first two possibilities however do not seem 
likely; one would require some as-yet unknown dependence on luminosity for the silicate 
feature, while the other would require an extraordinary level of serendipity to find nearly all 
of the sample in the same evolutionary stage. The third possibility is more plausible, but 
would require a 'top-heavy' IMF, for which there is currently no direct evidence. The fourth 
possibility is also plausble, and some level of selection bias is undoubtedly present in our 
sample as we are not selecting solely on the basis of IR luminosity, but selection bias likely 
plays only a minor role, for the reasons given above. We tentatively propose therefore that 
some combination of the final two possibilities is the most likely origin for most (but not all) 
of this effect. It is worth noting that these conclusions may also apply to sub-mm selected 
ULIRGs, given the similarity between their respective IRS spectra. 



^Such extended star formation would howev er still give rise to PAH emission, see for example the 'cirrus' 
models of Efstathiou &: Rowan-Robinson 20031 
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5. Conclusions 

We have presented mid-infrared spectra, taken using the IRS onboard Spitzer, of 32 
objects selected on (1) an optical band magnitude fainter than = 23, (2) excess flux in 
the IRAC 4.5/Ltm passband compared to the other three IRAC channels, and (3) a 24/im flux 
in excess of O.SmJy. Our conclusions are: 

1 - All of the sample lie within 0.8 < z < 2.3, with more than 90% lying in a much tighter 
redshift range, characterized by {z) = 1.71 ±0.15. The rest-frame IR luminosities are consis- 
tent with the range 10^^-^ — IO^^'^Lq. Our selection method is therefore a straightforward and 
efficient way of finding large samples of distant ULIRGs confined to narrow redshift ranges. 
It avoids all of the observational challenges inherent in finding distant ULIRGs via their 
far-IR emission, and produces samples with redshift distributions akin to those of optically 
selected samples such as Lyman-Break galaxies. 

2 - Nearly all of the sample show two or more strong PAH emission features in their mid-IR 
spectra, and systematically lack strong silicate absorption. The individual spectra closely 
resemble the IRS spectra of high-redshift sub-mm selected ULIRGs. The star formation 
rates computed from the PAH luminosities are all extremely high, of order > lOOOM0yr~^, 
sufficient to power all of the IR emission. Coupled with the lack of a power law continuum 
in the IRAC channels, we infer that star formation is likely to be the dominant power source 
in most of the sample, though it is possible that some also contain a luminous AGN. 

3 - Our sample occupy a different region on the PAH 6.2//,m EW vs silicate strength plane 
to most local ULIRGs, and there is marginal evidence for enhanced PAH 6.2/im/7.7//,m and 
6.2/im/11.2/im ratios. Furthermore, the composite mid-IR spectrum of our sample resembles 
that of local starburst galaxies with IR luminosities of 10^*^ — IO^^'^Lq, rather than that of 
local ULIRGs. Though selection bias hkely plays a role, we propose that the most likely 
reasons for this are that the star formation in distant ULIRGs is extended over scales of a few 
Kpc rather than the sub-Kpc starbursts seen in local ULIRGs, and/or occurs in unusually 
gas-rich environments. 
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Table 1. Observations summary 



ID Galaxy AOR key RA (J2000) Dec 3.6/^m 4.5AJm 5.8^lm 8^lm 24/im 
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Note. — All fluxes are given in fi]y. Errors on the IRAC fluxes are of order 5%, while the errors on the MIPS 24/im fluxes 
are of order 10%. A '— ' indicates an upper limit. Limiting fluxes are given in i|2.1l 



- 26 - 



Table 2. Redshifts & bump classifications 



ID Bump=' Zphot^ Zirs" 
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''Bump classification, based 
on the IRAC fluxes in Table [T] 

^Phot ometric redshifts taken 
from I Rowan- Robinson et all 
liool). The objects with a '?' 
are based on three photometric 
bands, and should be treated 
with caution. 

'^Spectroscopic redshifts de- 
rived from the IRS spectra. The 
three measurements with a '?' 
should be treated with caution. 
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Table 3. Fluxes and equivalent widths for major PAH features, silicate strengths, and star 

formation rates 



ID 
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Flux 
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EW 
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Flux EW 
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Flux EW 
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1 


3.88 ± 1.13 


1.01 ± 1.14 


5.67 ± 1.53 


0.53 ± 0.26 


2.68 ± 0.58 


0.63 ± 0.20 


0, 


,00 


± 


0.45 


1352 


13.28 


2 


<4.13 


<0.70 


<5.82 


0.43 ± 0.31 


2.53 ± 0.70 


0.48 ± 0.19 


0, 


,35 


± 


0.53 


<1850 


<13.46 


3 


<4.49 


<5.39 


8.62 ± 1.97 


0.77 ± 0.28 


3.99 ± 0.95 


1.44 ± 1.01 


1, 


,07 


± 


1.00 


2324 


13.49 


4 


8.82 ± 0.99 


0.64 ± 0.21 


17.79 ± 1.40 


0.61 ± 0.10 


5.03 ± 0.45 


0.51 ± 0.07 


0, 


,17 


± 


0.15 


2067 


13.53 


5 


2.47 ± 0.95 


2.66 ± 14.15 


5.27 ± 1.35 


0.58 ± 0.23 


1.55 ± 0.77 


0.37 ± 0.32 


0, 


,24 


± 


0.62 


1194 


13.24 


6 


3.15 ± 1.12 


0.52 ± 0.43 


9.22 ± 1.82 


0.80 ± 0.47 


2.39 ± 0.44 


0.50 ± 0.16 


0, 


,68 


± 


0.37 


1108 


13.16 


7 


<4.40 


<3.16 
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9 
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0, 
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± 
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1160 
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10 
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<0.64 


7.52 ± 1.85 
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1, 


,19 


± 
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<13.44 


11 


<4.03 


<0.83 


<6.0 
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<1.79 


<0.34 


0, 


,02 


± 


0.40 
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<13.38 


12 


3.46 ± 1.05 


1.02 ± 1.39 


4.48 ± 0.84 


0.38 ± 0.10 


<1.86 


<0.41 


0, 


,67 


± 


0.43 


<1458 


13.37 


13 


4.72 ± 1.38 


0.63 ± 0.59 


4.62 ± 1.53 


0.22 ± 0.08 


2.30 ± 0.42 


0.25 ± 0.06 


0, 


,20 


± 


0.26 


1160 


13.27 


14 


4.91 ± 0.96 


1.26 ± 1.52 


6.84 ± 1.51 
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0.66 ± 0.29 


0, 


,67 


± 


0.62 


1991 


13.50 


15 










<4.37 


<1.49 


1, 


,89 


± 


0.43== 






16 


<5.07 


<1.23 


2.80 ± 0.93 


0.22 ± 0.08 
















<13.72 


17 


1.53 ± 0.54 


0.63 ± 1.09 


3.94 ± 1.84 


0.57 ± 0.33 


2.33 ± 0.00 


0.59 ± 0.00 


0, 


,01 


± 


0.44 


1014 


12.94 


18 


3.96 ± 1.49 


0.98 ± 1.01 


4.25 ± 1.36 


0.26 ± 0.09 


2.89 ± 0.69 


0.53 ± 0.19 


0, 


,45 


± 


0.43 


1706 


13.38 


19 


4.94 ± 1.06 


0.95 ± 0.75 


9.65 ± 1.98 


0.59 ± 0.18 


2.99 ± 0.65 


0.44 ± 0.14 


0, 


,27 


± 


0.51 


1562 


13.38 


20 


<4.16 


<0.81 


<6.00 


<1.50 


4.23 ± 0.52 


0.93 ± 0.24 


0, 


,40 


± 


0.45 


<3008 


<13.74 


21 


4.35 ± 0.98 


0.56 ± 0.35 


7.64 ± 1.42 


0.54 ± 0.18 


1.45 ± 0.38 


0.24 ± 0.08 


0, 


,58 


± 


0.31 


1423 


13.42 


22 


5.85 ± 1.10 


7.34 ± 40.91 


8.28 ± 1.59 


0.58 ± 0.21 


3.89 ± 0.66 


0.71 ± 0.23 


0, 


,54 


± 


0.52 


2009 


13.48 


23 


<4.86 


<0.49 


4.68 ± 0.76 


0.51 ± 0.15 
















<13.77 


24 






8.16 ± 1.43 


0.31 ± 0.07 


8.83 ± 1.29 


3.32 ± 3.05 


0, 


,83 


± 


0.64 






25 


<3.10 


<0.28 


<5.00 


0.25 ± 0.11 


2.39 ± 0.68 


1.06 ± 0.66 


1, 


,34 


± 


1.34 


<1253 


<13.22 


26 


8.39 ± 1.21 


0.65 ± 0.22 


7.22 ± 1.47 


0.28 ± 0.07 


5.13 ± 0.75 


0.66 ± 0.15 


0, 


,51 


± 


0.35 


2424 


13.59 


27 


<3.31 


<0.65 


<4.80 


0.41 ± 0.22 


<1.30 


<0.26 


0, 


,01 


± 


0.43 


<1281 


<13.35 


28 


3.49 ± 1.19 


0.47 ± 0.45 


5.97 ± 1.23 


0.31 ± 0.09 


2.96 ± 0.40 


0.35 ± 0.06 


0, 


,00 


± 


0.24 


1030 


13.10 


29 


6.22 ± 1.10 


0.59 ± 0.25 


10.56 ± 1.57 


0.45 ± 0.09 


2.86 ± 0.55 


0.24 ± 0.05 


0, 


,73 


± 


0.27 


1929 


13.53 


30 


4.50 ± 1.12 


1.08 ± 1.06 


5.77 ± 1.62 


0.33 ± 0.11 


2.40 ± 0.50 


0.33 ± 0.09 


0, 


,37 


± 


0.32 


1533 


13.39 


31 


4.81 ± 0.93 


2.62 ± 5.95 


9.36 ± 1.18 


0.77 ± 0.31 


1.64 ± 0.49 


0.20 ± 0.07 


0, 


,85 


± 


0.47 


1370 


13.40 


32 


5.80 ± 0.90 


0.36 ± 0.10 


4.40 ± 0.96 


0.26 ± 0.07 


<1.85 


<0.19 


0, 


,02 


± 


0.50 


<2017 


13.60 


Note. — Fluxes are given in units of 10 ^^W cm ^ and equivalent widths are given in /xm. A '— ' indicates the feature lies outside 
the bandpass. 

"Star formation rates for those sources where both the 6.2/im and 11.2/jm PAH features lie within the IRS bandpass. The errors on 
all values are of order 50%, and assume that the star formation is similar in nature to that seen in local ULIRGs, but see i|4.3l 

''Rest-frame IR luminosity, computed using Equation [3] 



'^Measurement is inherently uncertain as the feature is not reliably identified, and only half in the bandpass 
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Fig. 1. — Individual spectra, plotted in the rest-frame using the spectroscopic redshifts in 
Table [2l The number in the top left of each panel corresponds to the ID number in the 
tables. The vertical dotted lines mark the wavelengths of the 6.2/im, 7.7/im, 8.8/im, 11.2/im 
and 12.7/im PAH features. 
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Fig. 3. — Left panel: Redshift histogram, using the redshifts derived from the IRS spec- 
tra. Nearly all of the sample is confined to a narrow redshift range of 1.4 < z < 1.9. 
Right panel: Comparison bet ween the IRS redshifts, and the photometric redshifts from 
Rowan-Robinson et al.l (120071 ). listed in Table [2l The dotted lines denote a deviation of 
0.06 in log(l-l-z), the boundary defined as a 'catastrophic failure' by iRowan- Robinson et al. 



( 120071 ) ■ Sources plotted as triangles have an uncertain IRS redshift. Sources plotted in grey 
have a photometric redshift derived from three bands. 
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Fig. 4. — Averaged spectra. The average of all 32 objects is plotted in black, with the grey 
shaded region indicating the la dispersion. The average of the ten objects with the largest 
PAH 6.2/xm EWs is plotted in red, and the average of the ten objects with the largest silicate 
strengths is plotted in blue. A 'P' denotes the wavelength of a PAH feature. 
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Fig. 5. — IRAC color-color plot. The 'wedge' is the selection region for sources that exhibit 
AGN-like IRAC colors, d efined empirically using IRAC observations of QSOs and Seyfert 
galaxies (ILacy et al.ll2004l ). 
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Fig. 6. — PAH luminosity v s PAH flux ratio plots. Our sample is plotted in black, and local 
ULIRGs ( ISpoon et al.ll2007l ) are plotted in green. The large orange point is the mean value 
for our sample. 
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Fig. 7. — PAH fl ux color-color plots. Our sample is plotted in black, local ULIRGs 



(ISpoon et al.l 120071) are plotted in green, and local 'classical' AGN (IWeedman et al.ll2005l ) 
and starbursts (iBrandl et al.ll2006l ) with IR luminosities in the range lO^'^L© - lO^^'^L© are 
plotted in red and blue, respectively. The large orange point is the mean value for our 
sample. 
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Fig. 8. — The 'Fork' diagram of ISpoon et al.l (120071 ). Our sources are plotte d in black 



l ocal ULIRGs are plo tted in green, while local starbursts (iBrandl et al.l 120061 ) and AGN 
(jWeedman et al.l 120051 ) with IR luminosities in the range 10^° — lO^^L© are plotted in blue 



and red respectively. The large orange point is the average for all the bump sources. Our 
sources occupy a small region, mostly the IC class. This contrasts with local ULIRGs, which 
are found in nearly all classes. 
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Fig. 9. — Comparison between our averaged rest-frame spectrum (black, with la errors 
denoted by the gr ey region) and averaged spectra for classes in the 'Fork' diagram of 
Spoon et al.l (120071 ). normalized at 14/xm; class lA (green), class IB (cyan), class IC (blue). 



class 2A (yellow), class 2B (magenta), class 2C (purple), class 3A(red), and class 3B (orange). 
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Fig. 10. — Comparison between our averaged spectra (black, \jir ~10 - 10 L©) 
and three different starburst galax y templates, all nor malized at 6/xm; the averaged spec- 

^Qio _ iQii-SL ^)^ Arp 

orange. 



t rum of the starburst g alaxies in iBrandl et al.l l2006l (blue, Lf b 
220 jArmus et al.l boOd green, L/r =1012.i4l^)^ ^nd M82 dsturm et al.l 
L/ij =10^°Lq, this spectrum was taken with the ISO-SWS and hence is of higher spectral 



2000 



resolution than the other spectra). Despite being approximately two orders of magnitude 
less luminous, the Brandl et al spectrum is an excellent match to our average spectrum. M82 
is also a reasonable match, except for slightly stronger apparent silicate absorption. Arp 220 
on the other hand is a very poor match, with much stronger absorption. 
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Fig. 11. — Spectral energy distributions of the eight objects in our sample with at least six 
optical-IR photometric data points (black circles) , plotted in the re st-frame. Overlaid are the 
best-fit templates derived from the iRowan-Robinson et al.l (120071 ) code. Redshifts are fixed 
at the spectroscopic values, and are given on the right hand side. In each fit, the leftmost 
dotted line is the starlight component, the righmost dotted line is an M82 template, and the 
solid black line is the sum of the two. For two objects (15 & 31) we show an alternative 
fit to the IR emission from an Arp 220 template (dashed line). The IRS spectra have been 
superimposed as red lines, but are not used in the fitting. 



